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ABSRACT 


Predictions of the sound pressure field in a wedge-shaped ocean overlying a fast 
bottom have been generated by computer-implemented calculations based on a method 
of images model and two parabolic equation models. Comparisons of the results show 
significant agreements and also some disagreements. The geometric shape of the wedge 
for the image calculation consists of a plane, sloping penetrable bottom and a plane, 
horizontal pressure release upper surface. This shape models a real continental shelf 

1 he shape used in the parabolic equation models is a pressure release plane surface 
above a conical bottom. Two PE models were studied. One. from a literature, is based 
on a wide-angle parabolic equation aiid providing a contour plot of TL as a function 
of position. The other, resident on the NPS computer, is based on an implicit finitc- 
difiorence (IED) algorithm. The results show that there is fair agreement among the 
difierent models. Comparison of the image method and the second PE model shows a 

2 to 3 dB dilference in transmission loss near the surface and good agreement deeper, 
'fhc transmission loss predicted by the first PE model diJfers from both other models by 
4 to 5 dB near the surface and at the middle depths. Near the bottom they all agree well. 
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I INTRODUCTION 


Acoustic propagation in an ocean with a sloping bottom is of 
growing concern both for theoretical reasons and because of the 
importance of knoAving the performance of acoustic sensors located 
over the continental slope. This problem has been studied by many 
scientists. In this research we compare and contrast the results of 
three different models, predicting the pressure field in a wedge- 
shaped ocean with a pressure release surface and an acoustically fast 
bottom. One image method and two PE models were selected. 

The computer model for the image method was written in Basic 
for use on a PC and in Fortran for use on the NPS IBM 3033 computer 
main frame. The majority of the computations for the image method 
was done using the Fortran version. 

A PE model contained in "Ocean Exec" library on the NPS IBM 
3033 computer was tested, but found to be too cumbersome and time- 
consuming. Modeling a sloping bottom with acceptable accuracy 
required 100 pairs of range-depth input which is tedious. In addition 
the program needs knowledge of certain parameters for a basement 
layer which is required to exist below the bottom. The output gives 
transmission loss at a fixed depth. To find the TL as a function of 
range for various depths between the bottom and the surface 
necessitated running the program many times which would have 
required an unacceptable amount of computer time. For these reasons, 
this PE model was not used for this work. 
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Two other PE models were found to be more acceptable. One was a 
model used at SACLANT by Tindle and Jensen iRef. 1] which we will 
refer to as PESAC. The other was a model based on an implicit finite- 
difference (IFD) algorithm implemented at NPS by Jaeger [Ref. 2]: we 
refer to this model as PEIFD in what follows. PESAC was not available 
for use, but some of its output was in the referenced article in the 
form of a contour plot of TL vs position in and below the wedge for a 
stated set of water and bottom parameters. The PEIFD, being available 
at NPS, was used to generate output for the same set of parameters. 
The results of the image method model and the two PE models are 
compared in the analysis that follows. 
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n. DEVELOPMENT 


A BIAGE THEORY 

In 1978. with the aid of a computer model based on the method of 
images, Coppens, Sanders, loannou, and Kawamura [Ref. 3] predicted 
and measured the pressure amplitude and phase in the upslope 
direction along the bottom of a wedge-shaped fluid layer overlying a 
fast fluid bottom. Baek [Ref. 4] used the same model in 1984 to predict 
pressure amplitude and phase everywhere within the wedge in the 
upslope direction. In the same year, LeSesne [Ref. 5] implemented a 
n odel developed by Coppens and Sanders which is not limited to up 
or downslope direction. In the method of images, boundaries are 
replaced by images of the source, placed geometrically to represent 
the propagation path between source and field point and given 
amplitudes and phases to reproduce the effects of interactions with 
the boundaries. For a wedge-shaped duct, images lie on a circle 
centered at the apex of the wedge. The source and each of the images 
radiate spherical waves of appropriate phases and amplitudes. The 
phase coherent summation of these waves yields the total pressure 
and phase at any field point in the wedge. In the computer program 
the number of images is made finite by insisting that the wedge angle 
P is an Integral submultiple of 180°. Then, the number of images 
(including the source) in upper or lower half-spaces is 

N = (180/p) (1) 
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Figure 1. Image structure for a wedge shaped duct (side view) 


The geometry of the image method is represented in Fig. 1. Figure 


2 contains the definitions and symbols for source and receiver 
parameters. 








surtace 



3 = wedge angle 

R] = scaled distance from apex to the source 

R 2 = scaled distance from apex to the receiver 

Yo = scaled distance along the shore between the 

receiver and the source 

Y = source angle measured upward from the interface 

5 = receiver angle measured upward from the interface 

Di = the ratio of medium density to the bottom density 

OC = the ratio of medium (water) sound speed to the 

bottom sound speed 

6c = critical angle for the bottom 

X = scaling distance 

a/k 2 = loss term in the bottom (used in the reflection 

coefficient) 


R(n) = range from the nth image to the receiver 
Figure 2. Definitions and notations of source and receiver parameters 
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For the wedge-shaped duct, sound will propagate from the source 
to the receiver through paths that can be divided into two families: (1) 
direct, surface reflected, surface then bottom reflected, surface- 
bottom-surface. and so on; (2) bottom reflected, bottom-surface 
reflected, bottom-surface-bottom, and so on For example, the bottom- 
surface-bottom path to the receiver in Fig. 3(a) can be interpreted as 
the straight line connecting the 3rd image with the receiver in Fig. 
3(b). 

The angle 6n of the n^h image measured from the bottom is shown 
in Fig. 1, 


= p(n-l) Y. 

for n odd 

(2.1) 

= pn - Y, 

for n even 

(2.2) 


The pressure from an image can be determined knowing the 
distance between the image and the receiver and the total reflection 
coefficient. The total reflection coefficient is the product of all 
reflection coefficients along the path from image to receiver. With 
pressure one meter from the source normalized to unity, the pressure 
at a field point from a specific image is found by dividing this total 
reflection coefficient by the distance. 
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source 


° receiver 



APEX 



(b) 

Figure 3. Visualization of (a) real path and (b) equivalent path in 
image method 


The distance between the receiver and the n^h image, as shown in 
Fig. 4 is 

R(n) = a/r? + R 2 - 2 R 1 R 2 cos(en-5) + Vg 

for the upper group of images, cind 

R'(n) = VR? + ^2 - 2 R 1 R 2 cosOn + 5) + Yg 
for the lower group of images. 










From Fig. 2, each image interacts with the bottom and surface a 
number of times; so that the pressure at the receiver from a specific 
image is proportional to the product of the relevant reflection 
coefficients. To find a reflection coefficient, the angle of incidence 
must be determined. The greater the index n of the image, the more 
reflections, and one must find an incident angle for each of these 
reflections. 
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, Risin(en-mp) + R2Sin(mp-5) 
sm(9n„)- 


Figure 5. The angle of incidence 0nm of the upper n^li image on the 
plane 
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for the lower images. [The notation is 6 for the upper images and 0' for 
the lower images.) 

The formula used in the program for calculating the reflection 
coefficient from a lossy bottom Is derived in the Appendix. The spatial 
phase term for each image is calculated from the distance. The 
pressure from an image is determined by multiplying the spatial phase 
factor with the total reflection coefficient and dividing it by the 
relevant distance. The summation of the complex pressure from every 
image is the complex pressure at the arbitrary receiver position. 

In the pressure calculation for upper images, it is not necessary to 
calculate reflection coefficient for the first and second image. These 
two images don't interact with any bottom plane. The reflection 
coefficient for the first image is 1 and for the second is -1. 
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The resulting complex pressure from the upper family of images is 


p„ = X ^expHkR(n)).(-l)l''/2| fjR,e^) 

n=iR(n) ni=l (5.1) 

OROnm) = 1. for n = 1.2 

and from the lower images is 

Pl= I :^exp{-JkR'(n)) (-l)'"/2I fiR(0'^) 

n=iKlnJ in=0 (5.2) 

In these expressions n indicates each image, and N is defined in 
Eqn. (1). The index m indicates each bottom and the number M is the 
total number of the bottom planes the sound meets during propagation 
to the receiver. It can be seen that 

M = [(n-l)/21 (6) 

The total complex pressure distribution P is the sum of the Pu and Pi; 

(7) 


P=Pu+Pi 





a PE MODEL 

As stated above, two parabolic equation (PE) models were 
compared with the image method. The first was the PESAC model 
from Jensen and Tindle (1987) which is based on the wide-angle 
parabolic equation. The pressure field in the model has cylindrical 
symmetry about a vertical axis through the source position. The 
sloping bottom therefor forms a conical hill for downslope propagation 
and a conical bowl for upslope propagation. A line source was used for 
the majority of cases. In one downslope case they used a simple, 
standard Gaussian PE source, with a -3 dB half-width of 35°, a wedge 
angle of 10° and a frequency of 5 Hz. The contour graph for this case 
will be compared to the results of the image method and a second PE 
model. 

The other PE model. PEIFD, was drawn from the thesis of Jaeger, 
This uses the implicit finite-difference (IFD) equation solution method 
to solve the PE. This method is unconditionally stable and has the 
capability to incorporate horizontal or sloping interfaces. The Fortran 
program in this thesis, like the PESAC model, specifies the receiver 
position with range and depth rather than range and angle. The input 
data used in this Fortran program were the same as the input data for 
Jensen and Tindle. 
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HL COMPUTATION OF IMAGE METHOD 

For the environmental parameters, the water has density and 
sound speed and the bottom has density p 2 , sound speed C 2 . and 
attenuation a/k 2 . 

To verify spherical spreading close to the point source, source and 
receiver were positioned as in Fig. 6. With the receiver in the same 
angular plane as the source, the receiver range from the source was 
varied. The results, Fig. 7. verify that the pressure near the receiver 
has a 1/r dependence. In Fig. 7(a), the source is at a scaled distance 
Ri = 2 and the multiplication of the range and the pressure is close to 
1, as it should be. The relation between scaled and real distance will 
be illustrated in chapter 4 with a specific case. The very first point at 
unit range is inaccurate because of calculational round off. 

Figure 7(b), where the source is at the scaled distance of Ri = 40, 
shows the same dependence within calculational inaccuracy. In this 
graph there were no results at receiver ranges of 1, 2, 3 m because of 
data overflow when source is this close to the receiver. This overflow 
is due to the value of the R(n). From Eqn. (3.1). R(l), the distance 
from the first image to receiver (direct path), is very close to zero. 
This near zero condition is a denominator for the pressure calculation 
of Eqn. (5), thus causing an overflow. 
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soil 




receiver 


Figure 6. The positioning of the source and the receiver 


The next three figures show dB loss vs position. For Fig. 8, the 
source is close to the apex and the pressure field is relatively simple, 
corresponding to the excitation of mostly the lowest mode. When the 
source range is relatively small this mode will dominate. The contour 
plots in Fig. 9 and 10 show that placing the source further from the 
apex excites the second mode. Placing the source even further from 
the apex will result in higher and more complex modes being excited. 
All the distances in Fig. 8. 9 and 10 are scaled distances. It is easy to 
show that the nondimensional cutoff distances of the first few modes 
are approximately 3, 5, 7, 11,.. 
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0 2 4 6 B 10 


RANGE FROM SOURCE (M) 

(a) source distance Ri = 2 


LO 



0 2 4 6 8 10 


RANGE FROM SOURCE (M) 

(b) source distance Ri = 40 

Figure 7. The range R(l) dependence of the pressure in the image 
method. Frequency 5 Hz, Dl = 0.5, CC = 5/6, (3 = 10°, y = 5 = 5° 
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310NV 

Figure 8. Excitation of mode 1. p = 10°. Ri =2 and the rest of 
parameters same as in Fig. 7 


16 


RANGE FROM SOURCE 














i I 

8 9 

310NV 

9. Mode 2 in a wedge, Ri = < 
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IV. DISCUSSION 


A. IMAGE METHOD AND PESAC 

The contour plot from Jensen and Tindle (1987) is shown in Fig. 
12(a). The input data for this contour are given by Table 1. The input 
data for the image method was chosen to be the same. 

The contour lines in Fig 12(a). which indicate dB loss, show that 
propagation is almost entirely in mode 1. In this figure, the frequency 
was 5.0 Hz, source depth (SD) was 50 m, water depth was 100 m and 
the maximum range was 10 km. 

To compare this with the image method, geometric parameters 
from the PE model were translated to parameters for the image 
method using Fig. 11. Source distance Ri. source angle y- and 
horizontal distance x for image method can be found from the wedge 
angle, water depth and source depth of the PE model. From simple 
geometric calculation, x = 567.13 m. y = 4.96°. Ri = 569.33 m. 


TABLE 1. PARAMETERS FOR MEDIA USED BY JENSEN AND TINDLE (1987) 



density(17mh 

sound speed 

attenuation 

water 

Pi = 1000 

Cl = 1500 m/s 


bottom 

p2 = 2000 

C2 = 1800 

aX2 = 1 dB 
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X 



source depth 50 m 
water depth 100 m 


Figure 11. Conversion of source distance and angle for image 
method from Jensen and Tindle (1987) 

To input ranges into the image method program, they should be 
converted to a scaled range. To convert Ri and R 2 into the scaled 
range, the following formula is applied. 


2 sin 9c tanp 

where X is the scaling range 

This scaled range is used in the image method program. The 
relation between seeded reinge and real range is 

R(scaled range) = R(real) / X 
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The critical angle (0c) is calculated to be 33.56°. With this critical 
angle, a wedge angle of 10°. and a frequency of 5.0 Hz. X = 769.48 m. 
This gives a scaled distance for Ri of 0.74. A scaled distance of R 2 = 
13.2 corresponds to the map range of 10.154 km. (The distance R 2 is 
measured from source to the receiver directly in the image method 
but horizontally in the PE model.) 

The attenuation in the PE model of aX 2 = 1 dB is equivalent to 
a/k 2 = 0.0183 for the image method since 

a/k2 = aX2/(8.7 x 2 tc) 

The transmission loss contours from the image method are shown 
in Fig. 12(b). These two contours show very similar shape. Figure 13 is 
an overplot of image-method results on the contour plot of the PE 
model. Near the surface and in the mid-depth level, the image method 
has 3 to 5 dB less loss, while they agree well near the bottom. 
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B. IMAGE METHOD AND PBIFD 

The PEIFD model is designed to generate transmission loss both in 
the bottom and in the water along horizontal range from the source at 
a given depth. The image method program was run in such a manner 
as to 5 deld comparisons between the image method and this PE model. 

The input file used for this purpose is shown in Table 2. The 
number in column 1. row 1 is the frequency in Hz and column 3, row 
1 is the receiver depth in m which was varied to probe the entire 
sound field, both in the water and in the bottom. Rows 3 and 4 define 
the sloping bottom and the definition of the remaining terms can be 
found in Ref. 2. 

In the image method, to find the pressure as a function of 
horizontal range at constant depth, the receiver angle must be 
changed, as the range is changed, as shown in Fig. 14. 







TABLE 2. INPUT DATA FILE FOR THE PE 


5 

50 

500 

0 

4000 


10000 

0.0 

0.0 

50 

1000 

25 

0 

100 





10000 

1763 





-1 

-1 





1763 

1.0 

0.0 




0 

1500 





1000 

1500 





1763 

1500 





8000 

2.0 

1.0 

1800 




6000 


The graphs shown in Fig. 15 through 22 indicate the results of the 
PEIFD model with a solid line and the image method with a dashed 
line. The PESAC model is annotated with the 'o' connected by dashes 
for selected graphs. The graphs are dB loss versus horizontal range in 
km at a fixed depth and all the source and receiver parameters are the 
same as used in the Jensen and Tindle model. The difference between 
the image method and PEIFD is approximately 2 to 3 dB in shallow 
water and the two lines converge as depth increases. The PEIFD 
results in the bottom also have a similar shape to the PESAC results. 
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Figure 15. dB loss .vs. range, depth = 50 m 
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Figure 16. dB loss .vs. range, depth = 100 m 
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Figure 17. dB loss .vs. range, dept±i = 200 m 
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Figure 18. dB loss .vs. range, depth = 300 m 
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Figure 19. dB loss .vs. range, depth = 500 m 


30 


RANGE FROM SOURCE (KM) 




I ^ , , I 

09 08 001 

(aa) ssoi NOiiV0Vdoyd 


Figure 20. dB loss .vs. range, depth = 800 m 
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Figure 21. dB loss .vs. range, depth = 1000 m 
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Figure 22. dB loss .vs. range, depth = 1500 m 


33 


RANGE FROM SOURCE (KM) 







V, CONCLUSIONS 

The geometric implementation of the ocean wedge and its 
boundaries by the image method and the PE models are different. To 
validate the use of these very different geometric models it is 
important to know how well they agree with each other for specific 
geometries. This thesis was an attempt to compare the image method 
and two PE models for propagation direction down slope and A\ith 
only one mode excited. 

Between the contour plot of the image method and PESAC model, 
there is a 3 to 5 dB difference at mid-depths. This difference mav be 
due to different normalization techniques or to the basic difference in 
geometry, but this hypothesis cannot be verified because the PESAC 
program was not available. The plot, however, shoves a very similar 
shape for the two models. 

On the other hand, when the image method is compared with the 
PEIFD model, the results show good agreement overall, with slightly 
less agreement at shallow depths. This would suggest that the 
different geometry may not lead to different TL's in this range of 
comparisons. Overall, however, since the image method compares 
quite favorably with both PE models, all three are apparently good 
tools for underwater sound prediction for downslope propagation with 
only one mode propagating. It is still an unresolved question how the 
various models will compare for upslope propagation, cross-slope 
propagation, and the source positioned to excite more than one mode. 
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APPENDIX 


The reflection coefficient is given by 


—sin 6-—sin 0* 

V - P| 

—sin0 +—sin0* 

Pi ^2 

where 0 is incident angle 

0t is transmission angle 

subscript 1 denotes water and 2 bottom 


and 



sin0t =^1-(C2 /ciPcos^0 

then 



—sin0^ =-)/(ci /C 2 )^ -cos^0 
C2 


With lossy bottom 


Sl 

C2 


Cl 




and 

^ (ci / C2 (1 + Ja / k2 f = (ci / C2 f + 2 J(ci / C2 a / k2 
C2 

makes 

-?i-sin0t = ^/(ci /C 2 )^ -cos^0 + 2j(ci /C 2 )^a/k 2 
C2 
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Let 


-^sinG^ = -^a-Jb = (a^ + 


where 


._, -ib 

a = (ci/c2p-cos20 b = 2(ci/c2)2a/k2 a 


With trigonometric identity 


J<t>/2 _ 


1 


l + Vl + (b/a)^ 

2Vl + (b/a)2 ^ 


Vl + (b/a)2-l 


^l + {h/aY 




+ a^ +a 


■'1 


^^2 ..2 


+ a^ - a 


2-^b^+a^ 


Utilizing this 


V2-^sin0t = VVb^Ta^Ta - jVVb 

C2 


2 2 
+a^ -a 


—sin 0 --|= VVPTa^T^ + j-)= VVb^+a^ 


V = 


_ Pi 


V2 


V2 


—sin0 + -^VVb^ + ^ + a - j-i'\/Vb^ +a^ - a 

Pj V2 

This form of the reflection coefficient is used to find pressure 
field in the program code for the image method. 
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